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Gas-Surface  Vibrational  Energy  Transfer  in  the  Transient  Region 
of  a  Low  Pressure  Unimolecular  Reaction 

+  * 

B.  D.  Barton,  D.  F.  Kelley,  and  B.  S.  Rabinovitch 

Department  of  Chemistry  BG-10,  University  of  Washington 
Seattle,  Washington  98195 

Abstract 

A  technique  has  been  developed,  called  the  Variable  Encounter  Method, 
for  the  study  of  gas-wall  vibrational  energy  transfer  in  the  transient  region 
of  a  unimolecular  reaction.  The  method  was  applied  to  the  isomerization  of 
cyclopropane  to  propylene.  Molecules  at  room  temperature  were  introduced  into 
a  hot  reactor  and  the  accommodation  to  the  final  temperature  of  the  reactor 
was  monitored  by  measurement  of  the  amount  of  reaction  that  occurred  at 
various  intervals  during  the  energy  relaxation  process.  Specifically,  the 
amount  of  reaction  following  some  (variable)  mean  number  of  collisions,  m  , 
between  gas  and  reactor  wall  was  measured.  Temperatures  in  the  range  900K-1125K 
were  employed.  Reactors  having  m  values  of  2.6,  8.5  and  27.2  were  used.  The 
average  probability  of  reaction  per  collision,  Fc(m),  for  a  given  reactor, 
was  deduced  from  the  data  and  compared  with  a  theoretical  stochastic  calculation 
based  on  both  gaussian  and  exponential  models  for  the  energy  transfer  proba¬ 
bilities.  It  is  found  that  the  efficiency  of  a  seasoned  quartz  wall  is  greater 
than  gas-gas  collisions  of  substrate,  and  that  the  efficiency  declines  with 
Increase  of  temperature.  The  steady  state  Is  closely  approached  (90%)  in  a 
comparatively  small  number  of  gas-wall  collisions  -  10-20,  approximately. 
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Introduction 

Lindemann  inspired  the  modern,  systematic  study  of  unimolecular  reactions 
His  mechanism  proposed  that  molecules  become  vibrationally  excited  by  colli¬ 
sions  with  other  molecules,  thereby  establishing  a  population  of  critically 
energized  molecules  which  may  decompose  to  products  if  not  deactivated  by 

collisions.  From  this  simple  mechanism,  the  study  of  energy  transfer  in  uni- 

2 

molecular  reactions  evolved.  At  sufficiently  low  pressures,  the  rate  of  a 

thermal  unimolecular  reaction  becomes  the  rate  of  critical  energization  of 

reactant  molecules.  The  rate  of  energization  depends  upon  the  nature  of  the 
3 

collision  partners.  The  study  of  energy  transfer  is  important  both  for  prac¬ 
tical  and  theoretical  reasons.  In  thermal  systems,  such  study  has  almost  in¬ 
variably  been  carried  out  under  steady  state  conditions. 

The  present  report  describes  a  study  of  energy  transfer  in  a  unimolecular 
reaction  in  which  molecules  experience  collisions  only  with  the  wall  of  the 
reaction  vessel.  The  objective  was  to  study  the  energization  process  in  the 
transient  region  prior  to  the  establishment  of  a  steady-state  population  of 
energized  molecules.  A  stochastic  treatment  has  been  applied  to  unimolecular 

4 

reactions  by  Shuler  and  Rubin  and  by  Montroll  and  Shuler.  Treatments  of  the 

"mean  first  passage  time"  (the  average  time  required  for  a  random  walker  to 

reach  an  absorbing  barrier  along  an  energy  coordinate)  have  been  developed  by 
5  6 

Kim  and  Widom.  In  principle.  If  one  knows  the  matrix  of  transition  proba¬ 
bilities  among  molecular  energy  levels  and  the  initial  population  vectors, 
then  one  may  calculate  the  mean  first  passage  time  and  the  transient  behavior. 

There  has  been  very  little  experimental  work  done  in  the  transient  region 
Dove,  Nip,  and  Teitelbaun/  found  an  "incubation  period,"  corresponding  to 
ca.  3000  collisions.  In  the  shock-heated  N2O  decomposition.  In  an  early  VLPP 

O 

paper,  Benson  and  Spokes  claimed  that  accommodation  to  an  energy  of  60  kcal 


3 

mole"^  occurred  in  *»»  80  collisions,  but  the  basis  of  their  treatment  was 

inadequate  to  support  their  conclusions.  In  this  paper  we  describe  a  new 

method,  called  £he  Variable  Encounter  Method,  for  the  study  of  the  transient 

region,  and  its  application  to  cyclopropane  isomerization.  A  preliminary 

g 

account  of  the  method  has  been  given. 


Description  of  the  Method 


The  Variable  Encounter  Method  (VEM)  depends  on  the  random  entry  of  mole¬ 
cules  into  a  hot  reactor  from  a  cold  reservoir  flask,  under  mc/lecular  flow 
conditions.  A  molecule  leaving  the  hot  zone  experiences  a  large  number  of 
collisions  with  the  cold  wall  of  the  reservoir  before  re-entering  the  reactor. 

The  molecules  thus  automatically  recycle  themselves  through  the  reactor.  The 
mean  number  of  hot  collisions  that  a  molecule  makes  per  encounter  with  the  hot 
reactor,  m  ,  is  a  function  of  the  shape  of  the  hot  zone.  One  may 

easily  measure  rather  small  reaction  probabilities  per  encounter  by  using  large 
experimental  times. 

Certain  design  features  of  the  experiment  are  important.  The  ratio  of 
cold  surface  area  to  hot  surface  must  be  sufficiently  large  to  ensure  that  a 
molecule  is  cooled  between  encounters.  The  mean  free  path  must  be  much  larger 
than  the  reactor  dimensions.  The  pressure  must  be  low  enough  so  that  the  system 
is  well  into  the  fall-off  region  (ideally,  in  the  second  order  region).  The  mean 
free  flight  time  between  wall  collisions  in  the  hot  zone  should  also  be  as  large 
as  possible  compared  to  the  mean  reaction  time  of  the  activated  molecules.  Of 
course,  the  final  exit  flight  involves  a  longer  path  length  into  the  cold  flask. 
The  total  amount  cf  gas  and  the  percent  reaction  are  to  be  commensurate  with 
the  sensitivity  of  the  analysis  method. 

The  relation  between  amount  of  reaction  and  the  experimental  time  is 

-!n(Rt/R0)  =  kobsl  (1) 

where  Rt  is  the  amount  of  reactant  at  time  t,  and  RQ  is  the  initial  amount 
of  reactant.  The  experimental  rate  constant  k  is  related  to  the  average 
probability  of  reaction  per  collision,  Pc(m),  for  a  given  (m)  reactor,  and  to 
the  rate  of  entry  of  molecules  from  the  reservoir  into  the  reactor;  the  latter 
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5 

involves  only  simple  mechanical  and  gas  kinetic  quantities. 

The  mean  number  of  collisions  m,  per  encounter,  was  calculated  by  a  ] 

Monte  Carlo  method  modeled  on  that  of  Davis, ^  and  the  principle  of  the  calcula¬ 
tions  is  reviewed  in  ref.  11.  The  cosine  law  of  reflection  was  employed  and 
is  believed  to  be  the  appropriate  model  under  the  present  condition.  The 
cosine  law  distribution  is  characteristic  of  diffuse  as  opposed  to  specular 
reflections  at  a  surface.  The  seasoned  quartz  surface  employed  is  assumed  to 

be  microscopically  rough.  Microscopic  roughness,  as  well  as  high  tempera- 

12 

tures,  usually  lead  to  diffuse  scattering.  In  any  case,  since  initial  entry 
into  the  hot  reactor  does  follow  a  cosine  law,  even  the  extreme  assumption  of 
specular  reflection  thereafter,  for  the  reactor  shapes  described  in  the  Experi¬ 
mental  section,  does  not  greatly  perturb  the  calculated  values  of  m  (1-2  units 
in  the  range  m  =  10  to  30). 

Along  with  m,  the  distribution  of  collisions  and  mean  distance  between 


collisions  were  also  calculated. 


Experimental 


The  apparatus  consisted  of  a  2-liter  fused  quartz  reservoir  sphere  on 
to  which  were  blown  several  cylindrical  closed  finger  reactors  of  diameter 
3.2  cm.  The  quartz  bulb  was  connected  to  the  system  by  a  ground  joint,  which 
permitted  each  of  the  reactors  to  be  rotated  into  position  for  heating.  In 
order  to  heat  a  reactor,  it  was  imbedded  in  a  stainless  steel  block  surrounded 
by  a  temperature-controlled  furnace.  Several  thermocouples  were  cemented  to 
the  heated  finger  so  that  the  temperature  along  the  reactor  could  be  monitored. 
Temperature  fluctuation  along  the  length  of  the  finger  amounted  to  less  than 
10°  over  most  of  the  length. 

The  length  of  a  particular  finger  determined  the  mean  number  and  distri¬ 
bution  of  collisions  that  a  molecule  made  on  entering  the  reactor.  Two  types 
of  reactor  apertures  were  used;  the  first  being  the  unconstricted  opening  to 
the  reactor,  and  the  second  being  a  blown  circular  constriction  at  the  entrance 
to  the  cylinder.  Reactors  having  m  values  of  2.6,  8.5  and  27.2  were  used. 

Some  data  were  collected  in  an  m  =  107  reactor,  but  were  more  fragmentary; 
they  are  also  believed  to  be  less  reliable,  because  the  reactor  was  apparently 
not  as  well  seasoned,  and  are  not  included  here. 

The  relative  sizes  of  the  quartz  sphere  and  the  aperture  of  an  attached 
reactor  were  such  that  the  internal  surface  area  of  the  sphere  was  more  than 
one  hundred  times  the  aperture  area.  Thus,  after  an  encounter,  molecules  were 
cooled  before  re-entering  the  hot  reactor. 

The  quartz  sphere  was  attached  to  a  glass  vacuum  system  that  could  be 
pumped  below  10"®  torr.  Provision  existed  for  gas  handling  and  measurement, 
sample  Introduction,  and  transfer  to  the  analytical  system. 

Reaction  mixtures  were  analyzed  by  gas  chromatography.  A  squalane  (28% 
on  Chromosorb  P)  and  phenyl  Isocyanate  -  Porasil  C  series  column  arrangement 
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was  used.  In  this  arrangement,  propene  was  eluted  before  the  parent  cyclo¬ 
propane,  with  a  concomitant  advantage  in  sensitivity. 

Cyclopropane  (Matheson  reagent,  99.0%  min.)  was  used,  usually  without 
further  treatment  for  most  experiments.  Analysis  revealed  0.3%  propene  as 
an  impurity.  For  some  experiments  involving  low  percent  reaction,  the  pro¬ 
pene  was  reduced  below  0.05%  by  passage  through  a  HgSO^-H^SO^  column  (1  cc 
of  a  20%  HgSO^-20%  H2S04  aqueous  solution  per  gram  of  Chromsorb  W). 

To  perform  a  run,  a  measured  aliquot  of  cyclopropane  was  admitted  into 
the  2-liter  quartz  bulb  for  a  given  time,  after  which  the  reaction  mixture  was 
pumped  out  through  a  liquid  nitrogen-cooled  trap.  The  time  required  to  pump 

out  the  vessel  to  1/e  of  the  initial  pressure  (~  10  sec)  was  entered  into  the 

-4 

run  time  with  suitable  weighting.  Run  pressures  in  the  bulb  ranged  from  10 

_3 

to  10  torr  Hg.  Run  temperatures  ranged  from  900  K  to  1126  K,  and  reaction 
times  from  25  seconds  to  4  hours. 


***** 
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Results  and  Discussion 

Experimental  data  were  obtained  by  measuring  the  percent  isomerization 
as  a  function  of  time,  for  each  reactor  at  each  temperature.  In  most  cases, 
three  separate  runs  were  made  at  a  given  run  time,  and  three  or  four  run  times 
were  used  so  that  the  determination  of  the  experimental  rate  constant,  and  of 
the  value  of  Pc (m) ,  at  a  given  temperature  in  a  given  reactor  was  the  result 
of  9  to  12  runs.  Plots  of  the  log  of  the  fraction  of  reactant  remaining  vs  run 
time  gave  straight  lines.  An  illustrative  plot  is  shown  in  Fig  1  for  the 
reactor  m  =  8.5,  T  =  1012  K.  The  plots  were  found  to  be  straight  lines,  as 
given  by  eq.  (1 ) . 

The  primary  quantity  of  interest  in  the  present  study  is  the  mean  proba¬ 
bility  of  reaction  per  collision  In  a  given  reactor,  Pc(m)  (Table  I). 

It  is  a  quantity  that  is  closely  related  conceptually  to  stochastic  theories 
of  chemical  kinetics.  Plots  of  ^(m)  vs  temperature  are  shown  in  Fig  2,  for 
mean  collision  numbers  m  of  2.6,  8.5  and  27.2.  It  is  noted  that  the  values 
of  Pc (m)  tend  to  converge  at  lower  temperatures,  i.e.,  the  wall  becomes  a 
stronger  collider  at  lower  temperatures. 

The  corresponding  Arrhenius  plots  give  slightly  curved  lines;  there  is 
no  a  priori  reason  why  the  plots  should  be  straight  lines  in  a  non-equilibrium 
system.  Average  Arrhenius  values  are  given  in  Table II. 

The  gas-wall  energy  transfer  process  may  be  described  in  terms  of  a  proba¬ 
bility  matrix,  where  p..  is  the  element  for  a  molecule  going  from  the  to 

*  J  — 

energy  level  upon  collision.  Once  the  functional  form  (model)  of  the  down- 
steps  of  this  matrix  is  specified,  the  corresponding  upsteps  may  be  calculated 

in 

by  detailed  balance.  The  details  of  the  construction  of  these  matricles  sre 
given  In  ref.  13.  Two  different  models  for  the  probability  of  a  down step  of  size 
AE  were  used. 
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exponential :  =  A  exp(-AE/<AE>) ;  0  s  AE  <  9000  cm" 

=0  ;  AE  >  9000  cm"^ 

gaussian:  =  A*  exp[-(AE-AEmp)^/2o^];  0  s  AE  s  9000  cm"^ 

=0  ;  AE  >  9000  cm~^ 

where  A  and  A'  are  normalization  constants  which  also  include  the  upstep 

probabilities;  <AE>,  AEmp  (mp  signifies  most  probable),  and  o  are  parameters 

of  the  model.  The  limitation  of  9000  cm~^  is  a  practical  computational  feature 

invoked  to  make  the  matrix  of  more  tractable  dimensions.  The  details  of  how 

calculated  Fc(m)  curves  are  obtained  from  these  models  of  energy  transfer  are 

given  in  ref.  14.  Because  of  the  truncation  feature,  we  also  define  an 

effective  value  of  the  average  energy  downstep,  called  <AE'>. 

In  all  cases,  the  values  of  <AE>  or  AEm„  were  selected  to  fit  the  m  =  2.6 

mp 

reactor  data,  and  F  (m)  values  were  then  calculated  for  the  other,  larger,  m 

V# 

reactors.  The  reason  for  this  procedure  is  that  the  calculated  value  of  Pc(m) 

is  somewhat  more  sensitive  to  the  value  of  <AE>,  or  AE^,  for  smaller  m  values. 

The  values  of  <AE*>  and  AEmp  used  at  various  temperatures  are  given  in  Table  III. 

It  is  seen  that  the  amount  of  energy  transferred  upon  collision  is  quite  large 

(approaching  strong  collider)  at  lower  temperature  and  decreases  at  higher 

temperatures  (Table  III).  A  decrease  in  efficiency  with  temperature  has  pre- 

15 

viously  been  observed  in  neat  gas  experiments  where  the  AE  values  are  signi- 

9  14 

ficantly  smaller.  This  same  decrease  has  been  observed  in  other  VEM  work  * 
and  in  some  low  pressure  pyrolysis  studies.^®  The  AE  values  reported  here  are 
somewhat  higher  than  those  found  for  cyclopropane-d2  (Table  III).  The  differ¬ 
ence  Is  systematic,  but  within  experimental  error  at  the  higher  temperatures. 

The  experimental  procedures  used  here  In  early  work,  particularly  seasoning, 
may  have  been  somewhat  less  reliable  than  in  the  other  studies.  The  present 
value  is  noticeably  higher  at  900K  although  still  qualitatively  and  semi- 


quantitatively  consistent.  It  is  seen  from  the  lack  of  fit  at  900K  that  an 

exponential  model  is  not  as  appropriate  as  a  gaussian  model  for  very  strong 

1 3 

collisions.  This  is  a  well-founded  conclusion. 


Conclusions 


The  results  of  the  VEM  experiments  for  cyclopropane  indicate: 

a)  below  900  K,  vibrational  energy  accommodation  is  substantially  complete 

in  'v  10  collisions;  at  higher  temperatures,  larger  collision  numbers  are 

necessary  for  complete  accommodation;  b)  wall  collisions  are  more  efficient 
15 

than  gas  phase  collisions,  with  wall  collisions  approaching  strong  collider 

behavior  at  lower  temperatures;  c)  collisional  efficiency  and  <AE>  decline 

with  increasing  temperature,  supporting  similar  findings  from  the  homogeneous 

15 

cyclopropane-d2  system. 
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Table  I 


Average  Probabilities  of  Reaction  per  Collision,  PCW 


m 

TOO 

Pc(m){xl08) 

2.6 

907 

2.71 

993 

9.40 

1081 

40.4 

8.5 

891 

7.1 

895 

6.6 

931 

15.8 

1012 

65 

1117 

512 

27.2 

906 

11.0 

989 

86 

1126 

1650 

13 


Table  II. 

m 

2.6 

8.5 

27.2 


Observed  Activation  Energies 


E.fkcal  mole' 

a 


30 

39 

47 
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Table  III.  Average  Values  of  <AE’>  (cm'^)  for  Cyclopropane- 

Wall  Collisions 


T(K) 

900 

1000 

1075 

Gaussian 

4900 

3275 

3000 

(a=0.35  AE  3 ) 
mp 

320014 

300014 

277514 

Exponential 

b 

3100 

2850 

303514 

277514 

242514 

a)  AE  is  same  as  <AE'>  within  a  few  cm"  ^ . 

mp 

b)  Exponential  model  too  inefficient;  cannot  give  fit. 
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Figure  Captions 


Figure  1  Experimental  decay  curve  at  m  5  8.5,  T=1012K 

Figure  2  Dependence  of  Fc(m)  on  T,  for  reactors  of  varying  m; 
dotted  curves  are  simple  exponential,  and  dashed 
curves  are  simple  gaussian  calculated  curves  for 
values  of  <AE>  and  AEmp,  as  in  Table  III,  fitted 
to  the  m  *  2.6  curve. 


0.0 


